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Abstract 



cr 

Vh^ We consider dynamics of a flat anisotropic multidimensional cosmological model in Gauss-Bonnet 

gravity in the presence of a homogeneous magnetic held. In particular, we hnd conditions under which 
the known power-law vacuum solution can be an attractor for the case with non-zero magnetic held. We 
^ | also describe a particular class of numerical solution in (5 + l)-dimensional case which does not approach 

| the power-law regime. 

q 

c4 1 I. Introduction. 

(N ■ 

The fact that structure of initial cosmological singularity can be rather complicated have been recognized 
since 60-th of the last century when the conception of BKL chaos have been presented [1]. It appears that 
Kasner solution being a general solution for a vacuum Bianchi I Universe becomes unstable in the case 
of metric with positive spatial curvature (belonging to Bianchi IX class) and is replaced by a complicated 
sequence of transient "Kasner epochs". Later it was found that some classes of an anisotropic matter can 
induce similar type of cosmological behavior even in flat Bianchi I case. This can be shown for a magnetic 
field by LeBlanc [2] and for a general vector field by Kirillov [3]. The BKL analogs for the flat magnetic 
Universe is a typical behavior when off-diagonal terms in the metric in the frame determined by magnetic 
field are present. For a diagonal case see for example [4]. 

Possible generalisation of these results may involve studies in the framework of multidimensional cos- 
mology as well as alternative gravity theories. Analogs of magnetic field for multidimensional space-times 
have been studied in [13, 14, 15]. If we want to introduce some corrections to Einstein gravity, the first goal 
is to find a solution which replaces Kasner solution for a flat anisotropic Universe. In the present paper 
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we consider Lovelock gravity as a generalization of General Relativity. The principal feature of Lovelock 
gravity is that this theory keeps the order of equations of motion the same as in GR, while other theories 
(like popular now /(i?)-theory) results in increasing of the number of derivatives in equations of motion. 
Another important property of Lovelock gravity is that it gives corrections to GR only in higher-dimensional 
space-time, so it is natural to consider it in the framework of multidimensional cosmology. As the number of 
non-Einstein terms in Lovelock gravity is finite for any given dimensionality of space-time, it is possible to 
consider regimes when these terms are not only small corrections to Einstein gravity. In particular, it is rea- 
sonable to expect that the highest order Lovelock term, consisting of highest power of curvature invariants, 
dominates near a cosmological singularity. In the present paper we consider the second Lovelock term, which 
is the famous Gauss-Bonnet combination. It is the highest possible Lovelock term for (4 + 1) and (5 + 1) 
dimensional spacetimes. As cosmology in the (4 + 1) dimensional case has some pathological features, we 
restrict ourselves by the (N + 1) dimensional case with N ^ 5 in the present paper. 

In late 80-th some vacuum solutions for a flat anisotropic Universe in Gauss-Bonnet gravity for N = 4, 5 
have been found [5, 6]. They replace Kasner solution of GR. Later these solutions have been rediscovered 
and verified for N = 6, 7 in [7]; after that this solution have been generalized for all N and also to the 
general Lovelock gravity [8]. The main goal of this paper is to address question of its stability near the initial 
singularity. It is known that in the presence of an ordinary isotropic matter these solutions are stable to the 
past if the matter has the equation of state with w < 1/3, otherwise the solution tends to isotropic one [9, 10]. 
For more interesting regimes near a singularity it is necessary to consider either anisotropic matter or curved 
geometries. As introduction of spatial curvature usually leads to very cumbersome equations of motion, we 
have chosen a magnetic fields as a possible source of instability. In this paper we restrict ourselves by the 
diagonal case. 

It is also known that in contrast to Einstein gravity certain initial conditions lead to exponential solutions 
instead to those of power-law behavior [9, 11, 12], we do not consider them in the present paper. 

The structure of the paper is as follows: In Sec. II. we describe the metric and matter content of the 
model studied, in Sec. III. the model is investigated in the framework of Einstein relativity. We present 
results on stability of Kasner solution in the presence of magnetic field. The same model in Gauss-Bonnet 
gravity is studied in Sec. IV. with presentation made as parallel as possible to presentation in the preceding 
section in order to compare Einstein and Gauss-Bonnet cases. In Sec. V. we describe a particular regime 
existing in the zone of instability of power-law regimes in (5+ l)-dimensional Gauss-Bonnet gravity. Sec. VI. 
contains a brief summary of results obtained. 

II. Preliminaries. 

In what follows we use a reference system chosen so that 1 : 



900 — 1) 9kk — — 



e 2ak(t \ 9lj = 0, ijtj 




Here and after Greek indices run from to N and Latin indices from 1 to N. 
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In this section we looking for the energy- momentum tensor of a pure magnetic field in (N + l)-dimensional 
space-time. In the general case components of the energy-momentum tensor of the electromagnetic field are 

T » = h ( F ^ F ^ + \^ F ^F aP ) (2) 
where F a p is the Faraday tensor. The components of the Faraday tensor obey the following equations: 

(3) 



V M F^ = 0. 
dF = 0. 



or, in more detail: 



where = As it is known 



+ r^ x F Xu + T^F^ = , 

(dF) afS ^ = 0, a < (3 < 7 , 



(4) 



therefore we may rewrite system (4) in the form 



| + T^ x F Xv + T^F^ = , 

\ T) a F pl - V p F ai + D 7 F Q/3 = 0, a < (3 < 7 

The components of the Riemann connection have the form: 



(6) 



r^ = r* oi = at®, I^ = d,-(t)e**« (7) 

Its other components are equal to zero. Since we consider the case of a homogenous space, Faraday tensor 
depends on time t only; then, in view of (7) the system (6) takes the form: 

F°* + F 0i J2a k = 0, 

k (8) 

F tJ =0 
Its solutions are 

~ a k 

F 0m = <j) m e k , 4> m = const , 
Fij = ipij, ipij = const 

Hereinafter we will be interested in the case of a pure magnetic field; so that (fi m = and, as a consequence, 

F 0m = F 0rn = (10) 

In what follows we will deal with diagonal energy-momentum tensor; in view of (1),(2),(9) and (10) it implies 
that 

Tj} = 0, ii ± v F V1 F^ = 0, ii ± v ^ e ~ 2ai ^ik = (11) 

l 

Functions e~ 2cil are linearly independent, therefore from (11) it follows that 



^ik = 0, i,k,l = l,N, i^k (12) 
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Each ipij is multiplied by all that ip^i, which has one of the indices k, I coincident with one of the indices 
Let us fix pair (i, j); the number of combinations, in which ipij is found, equals to 2(N — 2). Let ipij 7^ 0; 

then the other 2(N — 2) quantities ipik (V'j'fe); which are multiplied by ipij, must be equals to zero; indeed, if 

V'ifc 7^ (k i,j), then ipijipik 7^ 0, but that contradict (12). 

Continuing the argument, we find the number of zero magnetic field components: 

(N — 2) + (N — 2 N_1 ) N — 1 (N — I) 2 

2(N - 2) + 2(N - 4) + . . . = 2- ^-L = ^ '— for an odd dimensions (13) 

2 2 2 

(N - 2) + (N - 2#) N N(N- 2) , x 

2(N - 2) + 2(N - 4) + . . . = 2- ^ — = — i '- for an even dimensions (14) 



2 2 2 

ex como olpmanf c ami ale \r\ _ _ 



We took into account that the number of pairs of indices without the same elements equals to for an 
odd dimensions and to ^ for an even dimensions. 

The total number of components of the magnetic field is Af ^~ 1 ^ ; then the number \ °f non-zero magnetic 
field components is 

N(N-l) (N-l) 2 N-l r 
X = = — - — for an odd dimensions, (15) 

N(N-l) N(N-2) N „ 
X = = — for an even dimensions (16) 

In what follows they are exactly tpi2, V'34) ■ • • > V^n-i, 2n> n = 1> X that we set to be non-zero; other components 
of the Faraday tensor are assumed to be zero. Thus, components of the energy-momentum tensor of the 
electromagnetic field take the form: 

T„ = 0, M + v, (17) 

^0° = £ ^-i,*e- 2(0!K - 1+03i) , = -^L_ li2n e- 2 ^- 1+aa ») + £ ^-x.^e- 2 ^'-^^), n = L~x (18) 
i=i7x j=hx 

Hereinafter it will be convenient to use notations like this: 

<n>={2n-l,2n}, n = T^x (19) 

We assume that 

i ^< n > <=^> f ^ 2n - 1 A i ^ 2n, k =< n > <=^> fc = 2n - 1 V k = 2n (20) 

III. Stability of the Kasner solutions. 
III.l. Field equations. 



The action reads: 



S= lkl dN+lx V\det(g)\{C E + C m ) (21) 
C E = R, L m = F aP F aP (22) 
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where R and F a p are scalar curvature and Faraday tensor respectively. Here and after we use Planck units: 

c=l,G = -±-. = l, (23) 



N—2 ' 

mpi is the Planck mass. The gravitational equations is given by 

G$ = 8ttT£ (24) 

where 

G>t = K- U^R (25) 
are the components of the Einstein tensor. Taking into account (1) we get: 

G£ = 0, n + v (26) 
G% = . G n = Yl ( Si + ) + H n = ^ ( 2? ) 

As we consider the diagonal case here, then in view of (17), (18), (24), (27) for an even number of dimensions 
the gravitational equations can be written as 

Y, (a* + *?) + E = -^.-i.^-— 1+a9 »> + £ 4_ 1>2 ^- 2(a2 - 1+a2j) , n = (28) 

Y - E^-i.* e ~ a(a *- 1+aitf) = ( 29 ) 

i=Tx 

Note that each expression like (28) describes two equations with numbers 2n — 1 and 2n simultaneously 
(we use the notation (19) here). In the case of an odd number of dimensions there is one more equation in 
addition to these: 



Y (2* + 6?) + £ didj = £ 2] e- 2 ^ +a ^\ n=l, X (30) 

i,3^N J ' X 

The left side of the expression (29) is a first integral of (28). 



III. 2. Vacuum solution. 

When there is no any matter the equations (28)-(30) lead to: 



£(S 4 + h 2 ) + Y h i h j =0, n=l,N (31) 

Y OiOj = (32) 
i<j 

These equations has the following solutions [16]: 

a k = p k ln(t) + C k , C k = const, k = l,N, = 1, Y Pn = 1 ( 33 ) 

n n 

As a consequence, components of the metric tensor follows a power law: 

9kk(t) = e 2a ^ =C k t 2 P", C k = const, fc = MV (34) 
This is well-known multidimensional Kasner solution. 
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III. 3. Stability conditions. 

Metric ceases to obey a power law when there is a matter; but it turns out that the metric of the space filled 
with the magnetic field may get close to the Kasner metric and converge to it when moving towards the 
initial singularity. Namely, numerical calculations shows that there exists solutions of the equations (28) that 
get close to the Kasner solutions (33) and converge to them as time tends to the point t = 0. Considering 
an influence of the magnetic field as a perturbation we will search for solutions of the equations (28) in the 
form: 

a k (t) = a° k (t) + <p k (t), fc = L~iV (35) 

where a k is the Kasner solution (33), (f k G C 2 (K) is a perturbation. We will say that solutions a° k are 
asymptotically stable if 

lira a k (t) = a° k (t), k = hlV (36) 

t~>-0 

In other words, solutions a® are asymptotically stable if and only if 



lim ip k (t) = 0, k = l,N 



(37) 



Now we will find out conditions which specify asymptotically stable Kasner solutions. 

Proposition 1. Kasner solutions are asymptotically stable for t — > if and only if P2n-i+P2n < L n = 1,X- 
□ We consider the space of an even dimension; that of an odd dimension is treated similar way. 

1. Let {a\, . . . ,a^} be an asymptotically stable solution of the equations (31)-(32) given by (33); then there 
exists solution {ai, . . . , a^} of the equations (28), (29) and functions tpi, . . . , <pw S C 2 (R) such that 



a k (t) = a? k (t) + <p k (t), k=l,N 



limip k (t) = 0, k = l,N 



(38) 
(39) 



Substitution (33), (38) to the equations (28) and (29) leads to: 



j<k 
j,k^<n> 



j,k^<n> 



j,k^<n> 



-4>1 12 e ~ 2 (¥ , 2n-l+V2n)^2-2(p2n-l+P2n) _|_ ^ ^2 ^ e ~2(ifi2j-l+'fi2j)-j-2-2(p2j-l+P2j) 



J = 1,X 



It follows from (39) that 

lim(t<p k (t)) = 0, \im{t 2 (p k (t)) =0, k = h~N 
Let t ->• 0; in view of (33), (39), (42) equations (40)-(41) take the form: 



o = -ipl- 1 ,2j 2 - 2{P2n - i+P2n) + E ^-i.*' 



2-2(p 2j _ 1 +p 2j ) 



(40) 
(41) 

(42) 

(43) 
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0= £_Y>L 1>2 /~ 2(m - 1+P2<) (44) 

«=bx 

Equalities (43) and (44) are satisfied if and only if 

P2n-1 +P2n < 1, Tl=l,X (45) 

q.e.d. 

2. Let {a±, . . . ,a° N } be the Kasner solution given by (33) such that P2 n -i + P2n < L n = IjXj an d 
ifi(t), . . . ,(fN(t) £ C 2 (R) be a small deviations from that solution. We have: 

a k (t ) = a° k (t ) + <p k (t ) , k = MV (46) 

Let i = to > be the initial moment; we assume that 



\<Pi(to)\ « |a*(*o)|, \<Pi(to)\ < KM , i = l,N (47) 

Substitution (33), (46) to the equations (28) and (29) give us equations (40)-(41). In view of (47) we can 
neglect terms contained factors like ttpi, t 2 ipnpj in the lhs of the equations (40)-(41) and terms contained ipi 
in the rhs of that equations. Namely, let e be a small positive real number; then in the e-vicinity of the point 
to equations (40) can be written in the form: 

i 2 & = -^-^J 2 - 2(P2n ~ 1+P2n) + E ^f- 2(P23 - 1+P2j \ « = V* (48) 
or, in matrix notations, 

A$(i) = R(t), t £ (to — e, to + e) (49) 

where 

^| = 0, A) = 1, k ± j; $ fc (t) = fc (t), j, = TJlV (50) 

i?2n-l(t) = i?2n(t) = ^l-^J' 2 ^- 1 ^ ~ E 4-1, 2 ,r 2 ^^+^) , n = LX (51) 

It is easy to check that matrix A is nonsingular, so we obtain: 

$(*) = A _1 i2(t), (52) 

or, 

<Pk(t) = E ^kmt Xm , Vkm = const, A m = -2(p 2m -i + £>2m) , k = I, N (53) 

m=li 

Hence we deduce that 



E T ^ T t Am+1 , = E T^t Am+2 , fc = LiV (54) 



m=l,x rn=l,x 
It follows from P2m-l +P2m. < 1, = l,x that A m + 2 > 0; therefore 



lim = 0, k = l,N (55) 



and solutions (33) with P2m— l +P2m < 1> m = 1,X are asymptotically stable. 
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Examples. 

• (3+l)-dimensional space-time; in view of (15) we have a single non-zero component of the Faraday 
tensor; according to our convention (see Sec. II.) we denote this component by tpi2] under the Proposition 
1 it follows that stable Kasner solutions are characterized by the condition p\ + P2 < 1; it implies that 
P3 > — the space expands along the magnetic field. 

• (4+l)-dimensional space-time; in that case we have a couple non-zero components of the Faraday 
tensor: ip±2 and ^34; m accordance with the Proposition 1 it follows that stable Kasner solutions are 
specified by the criteria p\ + P2 < 1 and P3+P4 < 1. 



III. 4. Numerical calculations. 

The result presented above have been verified numerically. We have got a few hundreds numerical solutions 
with a random sets of initial conditions for each of the dimensions N = 5, 6, 7, 8 and found a lot of solutions 
that converge to the Kasner ones. Example of such solution is presented on the Fig. 1 (a). 

It should be note that looking for asymptotically stable solutions numerically using the synchronous time 
t has a number of drawbacks; particularly, numerical solutions of the equations (68)-(69) increase too rapidly 
impeding computations as t goes to zero. To avoid these problem we have introduced new time coordinate 
t by means of the relation 

-£«y(t) 

dr = e i dt (56) 

The change of variables from t to r results in transformation of logarithmic functions (33) into linear ones. 
Indeed, (33) and (56) leads to 

C , 
dr = —dt, C = e ' = const (57) 

Consequently, 

t = C ln(i) + C, C = const (58) 

a k (r) = ^r + C k - p k ^, k = TJV (59) 
Without lack of generality we can set 

d = . . . = C N = C = (60) 

Then, as was announced, 



a k (r)=p k r, k = l,N (61) 

Such a way, looking for asymptotically stable solutions numerically, we should expect to get quasi-linear 
solutions like ones that shown on the Fig. 1 (a). 
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Logarithm of components of the metric tensoi 




Determinant of M 



a) 



b) 



Figure 1: Numerical solutions for (5+l)-dimension space-time (Einstein gravity). The figure a) illustrates stability of the 
Kasner solution near the initial singularity. The figure b) shows decreasing of the determinant of the metric; it implies that any 
given element of volume dV = y] &et(g)\d N x tend to zero and therefore we go to the initial singularity. 

IV. Stability of power-law solutions in Gauss-Bonnet models. 
IV. 1. Field equations. 

Action of the theory under consideration is given by: 



1 



= ^ / d N+1 x^\tet(gj\{£ GB +C m ) 



C-GB — R a /3~fSR a ^ S — ^RapR aP + R Z i £-m — F a pF 



a/3 , P 2 



pa/3 



(62) 



(63) 



where R, R a p, R a ^5i Fap are the (N + l)-dimensional scalar curvature, Ricci tensor, Riemann tensor and 
Faraday tensor respectively. The gravitational equations has the form: 



H» = 8nT£ , 



(64) 



where 



Hi = 2RR£ - AR»RZ - AR^R" avp + 2R^R vaM - itfff (R aM5 R a ^ 5 - AR aP R a(3 + R 2 ) (65) 
are the components of the Gauss-Bonnet tensor. Taking into account (1) we get: 

HS = 0, ii^u (66) 

Hq = -12 ^2 hihjhkau H™ = -Ay^dj + °f) cij a k - 12 ^ didjd k d u n = l,N (67) 



i<j<k<l 



j,k^i,n 



i<j<k<l 
i,j,k,l^n 



For the diagonal case in view of (17), (18), (64), (67) for an even number of dimensions the gravitational 
equations can be written as: 



(a>i + a}) djd k + 3 djdj 



3<k 
j,k^i, <n> 



i<j<k<l 
i,j,k,l^<n> 
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2(a2j-i+a 2;) ) 



3 = liX 



n = l, X (68) 



£ didj-dfcd, + ^ £ ^.i,*^ 2 ^- 1 ^ = 0, (69) 

i<j<k<l i=Lx 

There is one more equation in addition to these in the case of an odd number of dimensions: 

^ai + af) £ Oja* + 3 £ Oi<W, = -i £ ^Li, 2l e" 2(a2l - 1+a2l) (70) 

i^iV j<k i<j<k<l i=~T~^ 

j,k^i,N i,j,k,l^N ,A 

The left side of the expression (69) is a first integral of (68). 
IV. 2. Vacuum solution. 

When there is no any matter the equations (68)-(70) lead to: 

^(oj + d 2 ) £ «jd fc + 3 Y djdjdfcdi = 0, n=l,N (71) 

ij^n j<k i<j<k<l 

£ aihja k ai = (72) 

i<j<k<l 

These equations has the following solutions [5, 6]: 

a k = p k ln(i) + C fe , C k = const, k = 1, N, £ PiPjPiPm = 0, = 3 (73) 

i<j<l<m n 

As a consequence, components of the metric tensor follows a power law: 

9kk(t) =e 2a *« = C k t 2p x, C k = const, fc = TTiV (74) 
Thus by analogy with (33)- (34) we will call solutions (73) and metric (74) Kasner-like for brevity. 

IV. 3. Stability conditions. 

As well as in the case of the Einstein gravity (see III. 3.), metric ceases to obey a power law when there 
is a matter; it appears however that the metric of the space filled with the magnetic field may get close 
to the Kasner-like metric (74) and converge to it when moving towards the initial singularity. Namely, 
numerical calculations shows that there exists solutions of the equations (68) that get close to the Kasner- 
like solutions (73) and converge to them as time goes to the point t = 0. We will looking for solutions of the 
equations (68) in the form: 



a k (t) = al(t) + <p k (t), k = l,N (75) 
where a? is the Kasner-like solution (73), ip k E C 2 (M). We will say that solutions a° are asymptotically 



stable if 



\ima k (t) = a k (t), k = l,N (76) 
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In other words, solutions a9 are asymptotically stable if and only if 



]hn<p k (t) = 0, k = l,N 

t — 5-0 



(77) 



Now we will find out conditions which specify asymptotically stable Kasner-like solutions. 

Proposition 2. Kasner-like solutions are asymptotically stable fort — > if and only ifp2n-i+P2n < 2, n = 1, X- 

□ We consider the space of an even dimension; that of an odd dimension is treated similar way. 

1. Let {a^, . . . ,a^r} be an asymptotically stable solution of the equations (71)-(72) given by (73); then there 
exists solution {ai, . . . , a^v} of the equations (68), (69) and functions tpi, . . . , (pry S C 2 (M) such that 



a k (t) = a u k (t) + <p k (t), k=l,N 



limip k (t) = 0, k = l,N 

t — 5-0 



(78) 
(79) 



Let us introduce the following notation: 



where ^4^ , . . . , ^.j^ , 0^ , . . . , Bi N are any indexed mathematical objects. Substitution (73), (78) to the equa- 
tions (68) and (69) leads to: 

Y PiiPi ~ 1) Y PjPk + 3 Y PiPjPkPl + 



i#<n> 



j<k 
j,kj^i 7 <n> 



i<j<k<l 



^2 Pi(Pi-l) Y P(jVk) + 2 Y P&i Y PjPk + 3 Y <P(iPjPkPl) 



i#<n> 



3<k 



3<k 



i<j<k<l 



+ 



j<fc 
J,fc^z, <n> 



j<fc 
J,fc^z, <n> 



+ 



X + X PjPk + 3 Y {^(iPkPl^j + ^(jPkPl) 



j<k 
j^k^i^<_n'> 



i<j<k<l 
i,j ,k ,i^<n> 



+ 



2 Y X + X yPi + Vi) Y P(j ( Pk)+^ Y Pii^WWl) 



3<k 
j,k^i,<Cn> 



j<k 
j,k^i,<n> 



i<j<k<l 
i,j,k,ly^<n> 



+ 



+t q 



Y + Y + 3 Y fiiVjvm 



3<k 



i<j<k<l 
i,j ,fc,Z^<n> 



,/,2 „-2(^2n-l+¥'2n)j.4— 2(p2n-l+P2n) 

^2n-l, 2n C 1 



2(V2j-l+¥'2j)^4-2(p2j-l+P2j) 



(81) 
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^2 PiPjPkPl+t ^2 ^(iPjPkPl) + t2 E (fiiPkPl)?] + ( Pi<P(jPkPl)) + 



i<j<k<l 



i<j<k<l 



i<j<k<l 
i,j,k,l^<.n'> 



i<j<k<l i<j<k<l 



It follows from (79) that 
In the limit i — >• in view of (73), (79), (83) equations (81)-(82) take the form: 



lim(^ fe (t)) = 0, lim(i 2 y3 fc (t)) =0, k = l,N 



(82) 



(83) 



= - lim 

4 t->o 



v-L-,^ 4 - 2 ^- 1 ^- EvW 



4-2(p2j_l+P2j) 







— lim -i/; 2 , , i 

12 i-> ^ ^2,-1,2, 



4-2(p 2l -l+P2i) 



i=l,X 



Equalities (84) and (85) are satisfied if and only if 

P2n-1 +P2n < 2, U = 1,X 



(84) 



(85) 



(86) 



q.e.d. 

2. Let {a^) • • • ; a %} be the Kasner-like solution given by (73) such that P2n-i + P2n < 2, n = l,x, and 
<Pi(t), ■ ■ ■ , <fN(t) G C 2 (1R) be a small deviations from that solution. We have: 



(87) 



afc(t) = ajfe(t) + ¥>fc(t), fc = l,JV 
Let t = to > be the initial moment; we assume that 

M*o)| < \a° k (t )\, \<fH(to)\ < |dg(* )| , i = MV 

Substitution (73), (87) to the equations (68) and (69) give us equations (81)-(82). In view of (88) we can 
neglect terms contained factors like tifi, t 2 ^pi^f>j, t^ipufijipk, t ■ l fi l Pi i Pk i Pl m the lhs of the equations (81)-(82) 
and terms contained cpi in the rhs of these equations. Namely, let e be a small positive real number; then in 
the e-vicinity of the point to equations (81) can be written in the form: 



f2 E & E p ^ pk = j 



or, in matrix notations, 



4-2(p 2 j-l+P2j) 



J = 1.X 



n = l, X (89) 



where 



A$(t) = i2(t), t e (t -e,t +e) (90) 



^2n-i = E Pip i> A 2n= E = n = l, X , k = l,N (91) 
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R2n-l{t) — R2n{t) — ~ 



2-2(p2j_l+P2j) 



3=1,X 



, n = l, X (92) 



Assuming that matrix A is nonsingular 2 , we obtain: 

&(t)=A- 1 R(t), (93) 

or, 

<Pk(t)= E S kmt lm , 8km = COIlst, 7 m = 2 - 2(p 2 m-l +P2m), k = 1 } N (94) 
m=l,x 

Hence we deduce that 

Mt)= E ^T^ +1 , <Pk®= E -^^ +2 , A; = lTiV (95) 

m=l,x 7n=l,x 

It follows from P2m-l + P2m < 2, m = 1, x that 7 m + 2 > 0; therefore 



lim = 0, fc = l,JV (96) 
t— 

and solutions (73) with P2m— l +P2m < 2, m = l,x are asymptotically stable. ■ 

Example. Let us consider (5+l)-dimensional space-time. According to (15) there are pair of non-zero 
components of the Faraday tensor; taking into account our convention (see II. for details) we denote them 
-012 and 034. Then under the Proposition 2 it follows that stable Kasner-like solutions are described by the 
conditions p\ + p2 < 2 and p% + P4 < 2. 

IV. 4. Numerical calculations. 

The result obtained have been confirmed by numerical calculations. As well as in the study of the stability of 
the Kasner solutions (see III. 3.) we have received a few hundreds solutions numerically with a random sets 
of initial conditions for each of the dimensions iV = 5, 6, 7, 8 and found numerous solutions that converge to 
the Kasner-like ones. Example of such solution is presented on the Fig. 2 (a). 

For the reasons mentioned above (see III. 4.) one should use special time coordinate for numerical calcu- 
lations. Let us introduce new time coordinate r by the following way: 

-s£«i(*) , . 

dr = e > dt (97) 

The change of variables from t to r results in transformation of logarithmic functions (73) into linear ones. 
Indeed, (73) and (97) give us: 

C -\X> G i 
dr = —dt, C = e = const (98) 



2 Matrix is singular iff its determinant is 0. Equation det(A(pi, . . . ,pjv)) = specifies a surface in (N — 2)-dimensional space 
of exponents pk that determine Kasner-like solutions (we have parameters pi, . . . ,pm connected by two relations: see (73)). 
Thus, excluding that surface from our consideration we lose no more than set of measure zero; as a consequence, our assumption 
has no significant effect on interesting result. 
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a) b) 

Figure 2: Numerical solutions for (5+l)-dimension space-time (Gauss-Bonnet gravity). The figure a) illustrates stability of 
the Kasner-like regime near the initial singularity. The figure b) shows decreasing of the determinant of the metric; it implies 
that any given element of volume dV = */] det{g)\d N x tend to zero and, consequently, we go to the initial singularity. 



Consequently, 

t = C ln(t) + C, C = const (99) 

a k (r) = ^T + C k - Pk ^, k = Tjf (100) 
Without loss of generality we can set 

d = . . . = C N = C = (101) 

Then, as was announced, 



a k {T)=p k T, k = l,N (102) 

Thus, searching for asymptotically stable solutions numerically we should expect to get quasi-linear solutions 
like ones that shown on the Fig. 2 (a). 



V. Oscillatory regime. 

Results of the preceding section pose an interesting question about the fate of trajectories in the instability 
zone. They can either reach stable Kasner-like attractor or represent some other type of dynamics different 
from Kasner-like behavior. In this section we show that the latter alternative can be realized. In our 
numerical calculations we have revealed a specific class of solutions of the equations (68). The discovered 
solutions possess the following properties: two of N functions a\, . . . oscillate and has opposite signs, 
the other ones vary slowly; owing to phase difference between the oscillating functions volume element 
det(g)\d N x oscillates too. In view of limited computational capabilities we have obtained numerical 
solutions of the equations (68) only for the cases N = 5,6, 7, 8; for each of that dimensions we have found 
out sufficiently large number of oscillatory solutions. Example of such solution is presented on the Fig. 3 (a); 
the corresponding oscillations of the volume element is shown on the Fig. 3 (b) . 

For N = 5 we have been able to derive an analytical approximation to such solutions on a certain 
interval of time. Below we describe the procedure of derivation of this approximation. Henceforth we use 
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Logarithm ol components of the metric tensor 



Logarithm of components of the metric tensor 





b) 

Figure 3: Numerical solutions for (5+l)-dimension space-time (Gauss-Bonnet gravity). The figure a) depicts oscillatory 
solution which does not tend to Kasner-like solutions both to the past and to the future. The figure b) shows the oscillatory 
behavior of the determinant of the metric; it means that any given volume element dV = y[det(g)|d x oscillates in this regime. 

time coordinate r introduced by (97). In the (5+l)-dimensional case equations (68) become: 



1 3 



4 [44 + 44 + 44] +444) = V* 5 J (tf 2 e-^- 2a * - 2^-^-^) 



i [44 + 44 + 44] + 444 

44 [4 + 4] + 44 [4 + 4 
44 [4 + 4] + 44 [4 + 4 



) = g e J=1 (V 1 



2 -2ai-2a2 n , 2 -2a3-2a4 
12 e z( / ; 34 e 



1 I E «i 



6 g ' 



= 1 lot? -2a 3 — 2a 4 



31' 



1 I E 



-e ' =1 (V>l 4 e~ 2a3 ~ 2a4 



2^ 2 2 e- 2ai - 2a2 ; 



2^ 2 2 e- 2ai - 2a2 ; 



(103) 

(104) 
(105) 
(106) 
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+ +aXK + 4])' = * % (i&e-**-** + ^e" 2 -- 2 -) (107) 

with the first integral 

a'^a'^ + ai444 + a'^a'^ + aiag^s + 02030405 = - * ' (ipl 2 e~ 2ai ~ 2a2 + ^| 4 e" 2a3_2a4 ) (108) 

The dash denotes derivative with respect to r. Analysis of numerical solutions shows that there exists To > 
(r = is accepted as the initial moment) such that 

o' 5 (r)«0, a^(T)^a' A (T), r > r (109) 

(see Fig. 6). In view of (109) equations (103)-(108) take the form: 



/ / I \' f ( I I /V / (ll l\' 9 ( I I /V 9 fnn\ 

0^304) = -, (030304) =-, (o 1 o 2 a 3 ) = -, [a^a^j =- (110) 



a i4[ a 3 + a 4] +a' 3 a' A [a' 1 + a' 2 ]\ = h (111) 



aia2 a 3 a 4 = 7 ( 112 ) 



where 



/ = ^ ie !^-§" _ 2K 2 et u -§ 1 ') , g = (-2 Kl e^ u + K 2 e^ v Y h = - ( Kl e* v -* n + K 2 e^ A (113) 

^ = 03 + 04, U = Ol + 02 (114) 

"1 = ^, «2 = ^f (H5) 
Equation (111) is the sum of the equations (110); equations (110) can be replaced by their equivalent: 

([oi+a'a]^)' = /, ([ai-c4]aX)' = 0, («i4[a' 3 + a 4 ])' = 0, (a'^ [a 3 - o 4 ])' = (116) 

In view of (109) the last equation from (116) is satisfied automatically, and the other equations become 

= 4/, (a'^v 1 )' = g, (K-o 2 ]^ 2 )' = (117) 



u v 



The last of equations (117) gives: 

C 
v 

so that 



a i ~ a 2 = if! C = const (118) 



Let us decompose functions /, g, h with respect to (|t> — |u) and (|u — |t>) up to the first order: 



2 4 
2k 2 ] - - [«i + 4k 2 ] « + - [«i + « 2 ] ^ (120) 



4 2 

« [« 2 - 2«i] +-[ki + k 2 ]u- - [« 2 + 4«i] u (121) 
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2 2 

h « — [kx + K2] + q [«i - 2k 2 ] u + - [«2 - 2«i] u (122) 

Numerical results show that such assumption is valid in a certain interval of time. Taking into account (119)- 
(122) we get: 

' u'v 



,,2^ -r „ n 8 r . . , . 16 



4 [«! - 2k 2 ] - - [k\ + 4k 2 ] u + — [kx + k 2 ] v (123) 



(C 2 1 \' 16 8 

(n') 2 -^ u'J =4[/e 2 -2/Ci] +y[Ki + K 2 ]u--[K 2 +4Ki]^ 



(124) 



Let us consider the equation (124). Numerical calculations show that in the regime under consideration 
v" /(v') A is the leading term in the lhs of the equation (124); then: 

C 2 v" 16 8 

= 4 [«2 - 2kx] + — [kx + k 2 ] u - - [k 2 + 4kx] v (125) 

We can obtain analytical solution of the equation (125) by omitting the term containing u in the rhs of (125): 

v" 8 

-jj = 4[«2-2«i] --[k 2 + 4kx]v (126) 
Solution of the equation (126) has the form 

v{t) = Tjy/r + A, rj, A = const (127) 
Let us consider the equation (123). We omit the term, containing v in the rhs of (123): 



u'v' 2 ' 



4[k! -2k 2 ] - -[kx + Ak 2 ]u (128) 



Under substitution of (127) into (128) we obtain: 

tu"(t) - u'(t) + i 2 r 2 u{r) = (t 2 (129) 
We have used the following notations: 

32[kx+4k 2 ] 16[«!-2« 2 ] 

€ = T~2 ' ^ = 2 

3r/ z 77^ 

Equation (129) has the following solution: 



u(t) = tCxJi \ ] + rC 2 Y2 | -i^— ] + ^, Ci, C 2 = const (131) 



3 \ 6 3 



£ 2 



Here J2 , Y2 are Bessel functions of the first and the second kind respectively. Going back to the initial 

3 3 

variables we get: 

C 

a[ - a' 2 = ^ > a 3 + «4 = u, «3 = 04, «5 = (132) 



Consequently, 



2C ~ ~ 

clx(t) — a 2 (r) = — ^r 2 + C, 03 (t) — C = aAr) = rjy/r + A, C,C = const (133) 
7/^ 



17 



and, finally: 



ai (r) -C t 2 + t 



02 (r) wC r 2 + r 



+ L>2 



where 



o 3 (r) « 7?^ + D 3 , 04 (r) « t/Vt + A, o 5 (r) « D, 
C ,Di,D 2 ,D 3 ,D = const 



(134) 

(135) 
(136) 

(137) 



To demonstrate oscillations of volume element -^/| det^Jld^x one should write down expression for the 
determinant of the metric: 



det{g)\ = e 2 (*i+-+«) = exp |4 {^yfF + r 



J D = J D + D 1 + D 2 + D 3 + A = const 



:i38) 



(139) 



Despite the fact that functions (134)-(136) are the solution of the simplified system (126), (128) they catch 
correctly qualitative features of solution of the full system (103)-(107) which can be shown by comparing 
functions (134)-(136) with numerical solution of the system (103)-(107). Figures 4,5 below depicts pairs 
(a,2i A2)', (04, A4) of numerical and analytical solutions; pairs (ai,A\); (03, A3) behave themselves analogously. 

This regime is not similar to Kasner-like one (Kasner-like solution has linear behavior on chosen time 
coordinate r (see IV. 4. for details), the leading asymptotic of the oscillatory solution has quadratic behaviour 
w.r.t. t). 
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Numerical and analitical solutions 





b) 

Fi gure 6i a) The figure shows that the derivative q!${t ) asymptotically tends to zero, b) The figure indicates that the difference 
a' 3 (r) — 014(7") asymptotically tends to zero. 
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VI. Conclusions 



In this paper we provided an analysis of stability of known power-law solutions in Gauss-Bonnet anisotropic 
cosmology in the presence of a homogeneous magnetic field in the case when metric is diagonal in the frame 
determined by magnetic field. The relative simplicity of equations of motion in this case allows us to find 
analytically a simple condition for the particular power-law solution to be stable. It appears that a set of 
initial conditions for which the power-law solution can not be an attractor near a cosmological singularity 
has non-zero measure (though does not coincide with the whole initial conditions space). The result has the 
same structure as in the Einstein gravity (particular sums of Kasner indices should be restricted by some 
number), the only difference is the exact value of this number. 

The fate of trajectories near singularity from this unstable domain requires future investigations. In the 
present paper we describe one particular class of trajectories found by numerical studies. The volume element 
experiences oscillations in the presented regime, making it qualitatively different from BKL oscillations. 

We can expect that some other regimes which do not tend to Kasner-like power-law regime exist in the 
system under investigation. Full description of them (the task which we leave for a future work) in this 
particular rather simple model can be useful in studies of cosmological singularity in Lovelock gravity in 
the presence of spatial curvature. The latter model usually lead to very cumbersome equations of motion, 
so understanding the nature of possible regimes in the instability zone for the Kasner-like solution in the 
present model would help in finding similar solution in more general and more technically complex situations 
of curved multidimensional geometries. 
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